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Abstract 

Thi 9 paper describes a joint NASA/Army re- 
search activity at the Langley Research Center to 
develop optimization procedures aimed at improving 
the rotor blade design process by integrating ap- 
propriate disciplines and accounting for important 
interactions among the disciplines. The activity 
Is being guided by a Steering Committee made up of 
key NASA and Array researchers and managers. The 
committee, which has been named IRASC (Integrated 
Rotorcraft Analysis Steering Committee), has 
defined two principal foci for the activity: a 

"white paper" which sets forth the goals and plans 
of the effort; and a rotor design project which 
will validate the basic constituents, as well as 
the overall design methodology for multidisci- 
plinary optimization. The paper describes the 
optimization formulation in terms of the objective 
function, design variables, and constraints. The 
analysis aspects are discussed, and an initial 
attempt at defining the interdisciplinary coupling 
Is summarized. At this writing, some significant 
progress has been made. Results are given in the 
paper which represent accomplishments in rotor 
aerodynamic performance optimization for minimum 
hover horsepower, rotor dynamic optimization for 
vibration reduction, rotor structural optimization 
for minimum weight, and integrated aerodynamic 
load/dynamics optimization for minimum vibration 
and weight. 


Introduction 

An emerging trend in the analytical design of 
aircraft is the integration of all appropriate 
disciplines In che design process (refs. 1 and 2). 
This means not only including limitations on the 
design from the various disciplines, but also 
defining and accounting for interactions so that 
the disciplines Influence design decisions simul- 
taneously rather than sequentially. Because the 
terms "integrated" and "discipline Integration" 
are frequently used Imprecisely, a definition of 
an integrated disciplinary design process is 
offered. Such a process is integrated if: 

(1) Information output from any discipline is 
expeditiously available to all other disciplines 
as required. 

(2) The effect of a design variable change pro- 
posed by one discipline on all other disciplines 
and the system as a whole is made known promptly. 


♦Senior Research Engineer, Member AIAA 
♦♦Chief, Aeromechanics Division, Member AHS 


Adhering to the above definitions is central to 
the plan to be described in this paper. The 
integrated approach has the potential to produce 
a be tter product as well as a better, more system- 
atic design process. In rotorcraft design (the 
rotor_in particular), the appropriate disciplines 
i nclude _ae rodynamic9 , dynamics, structures, and 
acoustics. The purpose of this paper is to 
describe a plan for developing the logic elements 
for hel icopter rotor design optimization which 
includes the above disciplines in an Integrated 
manner. 

Rotorcraft design is an Ideal application for 
integrated multidisciplinary optimization. There 
are strong Interactions among the four disciplines 
cited previously; Indeed, certain design parame- 
ters Influence all four disciplines. For example, 
rotor blade tip speed influences dynamics through 
the inertial and air loadings, structures by the 
ce ntrifug al loadings, acoustics by local Mach 
numb er an d air loadings, and aerodynamics through 
dynamic pressure and Mach number. All of these 
co nsider ations are accounted for in current design 
practice. However, the process is usually sequen- 
tial, not simultaneous, and often involves cor- 
recting a design late in the design schedule. 

Appli cations of rigorous and systematic 
analytical design procedures to rotorcraft have 
been increasing, especially in the past five 
years. Procedures have accounted for dynamics 
(refs. 3-8), aerodynamics (ref. 9), and structures 
(ref. 10). Generally, these applications have only 
cons i dered single-discipline requirements, al- 
though in reference 5, dynamic and structural 
requirements were considered together, and In 
reference 6, dynamics and aeroelastic stability 
we re combi ned. 

In early 1985, several occurrences led to an 
excellent opportunity at the NASA Langley Research 
Center to_ address the multidisciplinary design 
problem for rotorcraft. The Interdisciplinary 
Research Office was established and charged with 
the develo pment of integrated multidisciplinary 
optimization methods. Nearly concurrently, the 
Army Aerost ructures Directorate at Langley estab- 
lished the goal of improving rotorcraft design 
methodology by "discipline Integration." Close 
cooperation between the NASA and Army organiza- 
tions led to initial plans for a comprehensive, 
integrated analytical design capability. A group 
of NAS A/ Army researchers recently formed a com- 
mittee and began detailed planning for this activ- 
i ty* T he committee, designated IRASC (Integrated 
"Rotorcraft Analysis Steering Committee), has now 
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completed the bulk of the planning and has formu- 
lated the approach described in this paper. 


Sequence of Tasks 


Figure 2 depicts the general sequence of 
tasks that will lead to a fully integrated rotor 
blade aerodynamic/ dynamic/ structural optimization 
procedure which also accounts for acoustic and 
airframe dynamic influences. The dynamic optimi- 
zation work is building on the work described in 
references 5-7. The rotor aerodynamics activity 
has been separated into two parts. The first is 
aerodynamic performance optimization which is a 
continuation of the work described in reference 9. 
The second is an integration of aerodynamic loads 
analysis with dynamics - a procedure wherein the 
local airloads can be adjusted by varying the 
planform dimensions and twist of the blade to 
reduce dynamic response. A merger of the rotor 
performance optimization with the airload/dynamics 
optimization will yield a fully integrated 
aerodynamic/dynamic procedure. The rotor struc- 
tural optimization is a continuation of the work 
of reference 10, A merger of all the aforemen- 
tioned procedures, with the acoustic and airframe 
constraints interfaces, will lead to the fully 
integrated Phase 1 procedure. The resulting 
capability will be applied to a rotor test article 
to validate the procedures. 


The development of an integrated njltidlsci- 
plinary design methodology for rotorcraft is a 
three-phased approach. In phase 1, the disci- 
plines of blade dynamics, blade aerodynamics, and 
blade structures will be closely coupled, while 
acoustics and airframe dynamics will be decoupled 
from the first three but will be accounted for by 
effective constraints on the other disciplines. 

In phase 2, acoustics will be Integrated with the 
first three disciplines. Finally, in phase 3, 
airframe dynamics will be fully Integrated with 
the other four disciplines. In all three phases, 
systematically validated methods are the principal 
products of the research. 

This paper Is primarily concerned with the 
phase 1 activity; namely, the rigorous mathemat- 
ical optimization of a helicopter rotor system to 
minimize a combination of horsepower required at 
various flight conditions and hub shear transmit- 
ted from the rotor to the fuselage. The design 
will satisfy a set of design requirements includ- 
ing those on blade frequencies, autorotational 
inertia, aerodynamic performance, and blade struc- 
tural constraints. Additionally, the design is 
required to satisfy constraints imposed by re- 
sponse of the fuselage and also those constraints 
related to acoustics requirements* 


General Approach and Scope 
Development Strategy 

The general approach for the activity Is 
illustrated in figure I. In phase 1 the blade 
aerodynamic, dynamics, and structural analyses are 
coupled and driven by the optimizer. The optimi- 
zation of the blade aerodynamic geometry as well 
as the internal structure (spar, leading and 
trailing edge, ballast, etc.) is performed. The 
influences of the airframe dynamics and blade 
acoustics are accounted for in terms of design 
requirements (constraints) on the blade design. 
These requirements are described later in the pa- 
per. For a check on the efficacy of representing 
the acoustics requirements indirectly, the “final” 
design will be input to an acoustics analysis. 

The acoustics analysis calculates the acoustic re- 
sponse measures and derivatives of these measures 
with respect to the design variables. This infor- 
mation will be used to determine how well the 
design was able to satisfy the actual acoustics 
design requirements. 

The phase 2 procedure, wherein acoustics 
is fully Integrated with the blade aerodynamics, 
dynamics, and structural analysis, is also 
Illustrated in figure I. The design produced in 
phase 2 (when converged) will satisfy acoustics 
goals. Airframe dynamics In phase 2, as in phase 
l, is accounted for by effective constraints on 
the blade dynamics, aerodynamics, and structural 
behavior* Finally, In phase 3 airframe dynasties 
is integrated and the result is a fully integrated 
optimization strategy. 


Overall Problem Formulation 

This section of the paper consists of details 
of the integrated rotorcraft optimization problem. 
Included are descriptions of the following: the 

objective function (the quantity to be minimized 
for obtaining an optimum design); the design vari- 
ables (dimensions and other parameters of the de- 
sign); constraints (a set of behavioral or charac- 
teristic limitations required to assure acceptable 
and safe performance); and definitions of the 
Interactions among the disciplines. 


Objective Function 

The objective function will consist of a 
combination of the main rotor horsepower at five 
flight conditions plus a measure of vibratory 
shear transmitted from the rotor to the hub. 
Although several multiple objective function 
techniques are available (ref. li) one leading 
candidate is a linear combination whereby 

F - k 1 HP 1 + k 2 HP 2 + k 3 HP 3 

+ *4 HP 4 + k 5 RP 5 * k 6 S O) 

where F is the objective function 

kj through k^ are weighting factors 

HPj through HPe are required horsepower at 
various flight conditions 

S is the vertical hub shear 
A candidate set of flight conditions would be: 
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Flight 

condition 

Description 

Velocity 

(kts) 

Load 

factor 

1 

Hover 

0 

1.0 

2 

Cruise 

140 

1.0 

3 

High speed 

200 

1.0 

4 

Maneuver 

120 

3.5 

5 

Climb 

1000 fpm 
(VROC) 

— 


Blade Model and Design Variables 

Figure 3 Is a depiction of the rotor blade 
model to be used in the phase i optimization 
activity. Also shown in figure 3 are the design 
variables which are defined in table 1. The blade 
model may be tapered in both chord and depth. The 
depth is linearly tapered from root to tip. The 
chord is constant from the root to a spanwise lo- 
cation (referred to as the point of taper initia- 
tion) and may be linearly tapered thereafter to 
the tip* Design variables which characterize the 
overall geometry of the blade include the blade 
radius, point of taper initiation, taper ratios 
for chord and depth, the root chord, the blade 
depth at the root, the flap hinge offset, and the 
blade maximum twist. Tuning masses located along 
the blade span are characterized by the mass 
values and locations. Design variables which 
characterize the spar box beam cross section 
Include the wall thicknesses at each spanwise 
segment and the ply thickness at 0° and ±45®. 
Additional design variables include the number of 
rotor blades, the rotor angular speed, and the 
distribution of airfoils. 


Constraints 

As previously described, the phase 1 activity 
is based on integrating the blade aerodynamic, dy- 
namic, and structural analyses within the optimi- 
zation procedure. The acoustics and airframe dy- 
namics analyses are decoupled from the first three 
disciplines and their Influences are expressed in 
terms of constraints. Accordingly, the set of 
constraints Is made up of two subsets. The first 
subset consists of constraints which are evaluated 
directly from the first three disciplinary analy- 
ses and are a measure of the degree of acceptabil- 
ity of the aerodynamic, dynamic, and structural 
behavior. The second subset represents indirect 
measures of the satisfaction of constraints on the 
acoustics behavior and the requirement of avoiding 
excessive vibratory excitation of the airframe by 
the rotor. 

The constraints are summarized in table 2. 

The first two constraints are for aerodynamic per- 
formance and require that for all flight condi- 
tions, main rotor horsepower not exceed available 
horsepower and that airfoil section stall not 
occur at any azimuthal location. The next nine 
constraints address blade dynamics. The first 
requires that the blade natural frequencies be 
bounded to avoid approaching any multiples of 
rotor speed. The next five Impose upper limits on 
the blade vertical and inplane loads, transmitted 
hub shear, hub pitching, and rolling moments. The 
next three dynamic constraints are an upper limit 
on blade response amplitude, a lover limit on 
blade autorotat tonal inertia, and finally, the 
aeroelaatlc stability requirement. The structural 


constraints consist of upper limits on box beam 
stresses, blade static deflection, and blade twist 
deformation. The acoustic constraints are ex- 
pressed as an upper bound on the tip Mach number 
and an upper bound on the blade thickness to limit 
thickness noise; and an upper bound on the gradi- 
ent of the lift distribution to limit blade vortex 
Interaction (BVI) and loading noise. The airframe 
constraints are expressed first as a separation of 
the fundamental blade inplane natural frequency in 
the fixed system from the fundamental pitching and 
rolling frequency of the fuselage to avoid ground 
resonance; second as a bounding of the blade 
passage frequency to avoid the proximity to any 
fuselage frequency. 


Interdisciplinary Coupling 

Phase 1 of the effort will utilize several 
design variables which have historically been 
significant drivers of disciplinary phenomena. In 
addition, other variables are being included to 
provide other unexplored design opportunities. 
Table 3 shows an attempt to characterize the in- 
teractions among the disciplines through the 
design variables. For example, rotor tip speed 
has driven past rotor designs based solely on 
acoustics, performance, or dynamics. This vari- 
able also Influences blade structural integrity 
and fixed system response to transmitted loads. 
This provides the strong interdisciplinary cou- 
pling for tip speed shown in table 3. There are 
variables, such as blade twist, which can strongly 
Influence some disciplines, such as aerodynamics, 
while not perturbing others (e.g. , structures) and 
other variables such as a hinge offset which, 
heretofore, have not greatly Influenced conven- 
tional rotor design. 

A significant part of the current effort 
will not only explore the obvious strong design 
variable couplings, but will also address those 
variables which may provide design synergism for 
multidisciplinary design goals. This may provide 
a design key for missions which have not been 
accomplished with today’s rotorcraft. 


Implementation Method 
Organization of System 

The overall organization of the system to op- 
timize a blade design for aerodynamics, dynamics, 
and structural requirements is shown schematically 
in figure 4. In order to perform the aerodynamic, 
dynamic, and structural analyses indicated in the 
blocks in figure 4, it is first necessary to 
transform or "pre-process" the design variables 
Into quantities needed in the various analyses. 
Forrexaaple, the dynamic and structural analyses 
^oth need stiffnesses ~ El - and GJ, and laminate 
properties. The aerodynamic analysis needs lift 
and_drag coefficients for the airfoils used. 

The above Information is obtained by the design 
variable pre-processors which act as translators 
of the global design variables into local vari- 
ables needed In the analyses. The output of each 
analysis block, in general, serves two purposes. 
First, response-type output may be transmitted to 
another analysis block (e.g., airloads from aero- 
dynamics to dynamics); second, information is 
supplied to the objective function and constraints 
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block (e.g. , stress constraints from the struc- 
tural analysis). A key part of the procedure is 
the sensitivity analysis. This block corresponds 
to the calculation of derivatives of the con- 
straints and objective function with respect to 
the design variables- The derivatives quantify 
the effects of each design variable on the design 
and, thereby, identify the most important design 
changes to make enroute to the optimum design. 

The sensitivity information is passed to the 
optimizer along with the current values of the 
design variables, constraints, and objective 
function. The optimizer uses the information to 
generate a new set of design variables, and the 
entire procedure is repeated until a converged 
design is obtained. For our purposes, a design is 
converged when all constraints are satisfied and 
the objective function has reached a value which 
has not changed for a specified number of cycles. 


Optimization Algorithm 

The basic optimization algorithm to be used 
in this work is a combination of the general- 
purpose optimization program C0NM1N (ref. 12) and 
approximate analyses for computing the objective 
function and constraints. Because the optimiza- 
tion process requires many evaluations of the 
objective function and constraints before an 
optimum design is obtained, the process can be 
very expensive if complete analyses are made for 
each function evaluation. However, as Miura 
(ref. 3) points out, the optimization process 
primarily uses analysis results to move in the 
direction of the optimum design; therefore, a 
complete analysis needs to be made only occasion- 
ally during the design process and always at the 
end to check the final design. Thus, various 
approximation techniques can be used during the 
optimization to reduce costs. In the present 
work, the objective function and constraints will 
be approximated using piecewise linear analyses 
that consist of linear Taylor series expansions. 


Analysis Aspects 

The analytical tools must provide technical 
fidelity in phenomena predictions, as well as 
connectivity between disciplines. The areas 
of aerodynamics, dynamics, and structures will 
utilize codes to predict response, as well as 
sensitivity information. The constraint-providing 
disciplines of acoustics and airframe dynamics 
have the analysis task of defining the impact of 
the design on acoustic energy and fuselage 
response. 

The aerodynamic analysis for rotor perfor- 
mance prediction will include a hover momentum/ 
strip theory code for hover and climb applications 
(ref. 13). The CAMRAD program (ref. 14) will he 
used for forward flight and maneuver performance. 
In order to assure that the latest developments in 
inflow analyses are available, some modularity 
will be provided in the inflow modeling based on 
recent fidelity assessments. 

Rotor dynamics will utilize CAMRAD for forced 
response calculations. Finite element modeling 
(ref. 15) and the modified Galerkin technique in 
CAMRAD will form the tools for the dynamic tuning 


before the global analysis predicts the final 
blade loads, response, and rotor stability. 

The structural codes involve a combination of 
beam analysis and laminate analysis. The analysts 
(e.g., ref. 10) is applied to the blade planform 
model. The laminate analysis will be applied to 
one or more cross-section models. The beam model 
consists of equivalent stiffness and masses from 
which displacements and forces are computed. The 
Internal blade structure is represented by cross- 
section models to calculate resultant stresses 
associated with each beam model segment. The 
laminate analysis then uses these stresses to 
determine critical structure margins of safety. 

The effectiveness of imposing phase 1 
acoustic constraints will be quantified by using 
the WOPWOP code (ref. 16), with appropriate load- 
ing inputs from CAMRAD. Low frequency loading, 
thickness, and BVI noise will be generated from 
this analysis. 

Airframe dynamics constraints for phases 1 
and 2 will result from fixed-system frequency pre- 
dictions and will neglect hub motion. Phase 3 of 
the effort will involve finite element modeling 
and impedance tailoring to effect favorable rotor- 
body coupling in the design process. 


Validation Strategy 
Validation of Procedures 

The process of validating the optimization 
methodology involves substantially more than eval- 
uating the success of the final design. Specifi- 
cally, the analyses used in optimizing the rotor 
during Phase 1 will be examined for predictive 
fidelity and design technique validation. The 
usefulness of the basic tools involves not only 
accuracy of analysis, but also a reliable paramet- 
ric sensitivity capability. Several opportunities 
are currently available to assess the fidelity of 
the analyses. For example, rotor performance, 
dynamics, and acoustics predictions need accurate 
inflow distributions for various flight con- 
ditions. Recent experimental efforts (e.g., 
ref. 17) and code validations (ref. 18) are help- 
ing to provide confidence in the available inflow 
models. Rotor geometric design variable sensitiv- 
ity (e.g., effect of taper on performance), which 
was reasonably well-known for past rotor designs, 
is being re-examined in light of recent correla- 
tion anomalies for high-speed flight. Acoustic 
source mechanisms and modeling validity are also 
being examined (ref. 19), especially for paramet- 
ric sensitivity of the acoustic energy to rotor 
state. Structural coupling mechanics are being 
exploited in new rotor designs to assess the 
structural tailoring benefits while satisfying 
structural integrity requirements (ref. 20). 

Proof of the fidelity of design techniques is 
crucial to the overall design optimization effort. 
For example, aerodynamics and dynamics interact so 
strongly in rotor design that baaic aeroelastic 
tailoring efforts must be validated. Such a vali- 
dation effort is being undertaken at Langley, as 
well as other research centers (ref. 21). Also, 
because rotor speed is a strong driver for aero- 
elastic response, a program to assess variable RPM 
designs is underway at Langley. The objective of 
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this effort is to define the benefits and limita- 
tions of an aerodynamically and dynamically de- 
signed multi-speed rotor. In addition to design 
techniques which capitalize on the strong effects 
of certain design variables, small variances in 
other blade characteristics may Impede the practi- 
cal operation of even conventional designs. Hence, 
the ability to accurately predict even these sec- 
ondary phenomena Is Important for the design 
effort. For example, t rack-and-balance sensitiv- 
ity experiments and studies are being undertaken 
which can lead to practical design capability to 
eliminate blade-to-blade variability effects. 


Overall Design Validation 

For the overall phase 1 validation effort, 
the Langley team is defining a rotor task which 
requires maneuverability, speed, and efficiency 
(see table 4). Specifically, the rotor mission 
must be accomplished with minimum power and 
vibration while satisfying predefined acoustic, 
stability, and fuselage dynamics requirements. 
This validation activity is, in effect, a design 
project which will produce a rotor teat article. 

The assessment of the phase 1 design methods 
will involve model rotor hover and wind tunnel 
tests. The models (a baseline and an advanced 
design) will be aerodynamically and dynamically 
scaled. Provisions for varying key design param- 
eters are necessary to complete the validation 
process. In other words, the tests need to quan- 
tify not only the minima, but the gradients. 

The testing possibilities include a series 
of l/5-scale model rotors, mounted on a variable 
drive system and tested In hover and simulated 
forward flight in a tunnel which can eliminate 
many testing "excuses” such as inappropriate 
Reynolds, Mach, and Froude Numbers* The Langley 
TDT is the candidate facility for the major 
segments of the validation process. 


Results Obtained to Date 

Progress has been made in the areas of aero- 
dynamic performance optimization, optimum place- 
ment of tuning masses for vibration reduction, 
structural optimization, and integrated aerody- 
namic load/dynamic optimization. Selected results 
from these activities are highlighted In this 
portion of the paper. 


Results - Aerodynamic Performance Optimization 

A Mathematical Programming technique (ref. 9) 
has been developed to minimize the hover horse- 
power for a helicopter with a specified design 
gross weight operating at a specified altitude 
and temperature (fig. 5). A conventional design 
approach Is usually a two-step iterative method. 
The first step Is design for optimum hover per- 
formance by varying taper ratio, point of taper 
initiation, and twist until the rotor blade con- 
figuration with the lowest hover horsepower Is 
obtained. In the second step, this best hover 
design is modified by changing the root chord to 
meet forward flight and maneuverability require- 
ments. The Mathematical Programming approach used 
the same performance analyses as the conventional 
approach, but coupled a general-purpose optimiza- 


tion program to the analyses. The conventional 
and Mathematical Programming approaches have been 
used to define the blade configuration which pro- 
vides the lowest hover horsepower and satisfies 
forward flight and maneuverability requirements. 
Figure 5 summarizes results for the final design 
variable values and the main rotor horsepower 
required for hover from each approach. The 
Mathematical Programming approach produced a de- 
sign with more twist, a point of taper Initiation 
further outboard, and a smaller blade root chord 
than the conventional approach. The Mathematical 
Programming design required 25 less hover horse- 
power than the conventional design. Most signif- 
icantly, the Mathematical Programming approach 
obtained results more than ten times faster than 
the conventional approach. 


Results ~ Optimum Locations of Vibration 
Tuning Masses 

The objective of this work (described in 
ref. 22) is to develop a method for optimally lo- 
cating, as well as sizing, tuning masses to reduce 
vibration using formal mathematical optimization 
techniques. The design goal is to find the best 
combination of tuning masses and their locations 
to minimize blade root vertical shear without a 
large mass penalty. The method is to formulate 
and solve an optimization problem in which the 
tuning masses and their locations are design 
variables that minimize a combination of vertical 
shear and the added mass, with constraints on fre- 
quencies to avoid resonance. Figure 6 shows an 
arbitrary number of masses placed along the blade 
span. The optimization strategy reduces the 
oscillatory shear as a function of time during 
a revolution of the blade. 

Results have been obtained wherein the above 
strategy was applied to a rotor blade considering 
multiple blade mode/multiple harmonic airloads. 

The example problem is a beam representation of an 
articulated rotor blade. The beam Is 193 inches 
long with a hinged end condition and is modeled 
by 10 finite elements of equal length. The model 
contains both structural mass and lumped (non- 
structural) masses. Three lumped masses are to be 
placed along the length of the beam. The strategy 
was applied to a test case of two modes responding 
to three harmonics of airload. Figure 7 shows for 
the initial and final designs, the shear s plot- 
ted as a function of the time and azimuth for one 
complete revolution of the blade. The peaks on 
the initial curve have been reduced dramatically. 
For example, the maximum peak oscillatory shear 
for the initial design is 78.00 lbf , and for the 
final design, the maximum peak is 0,60 lbf. 


Results - Rotor Structural Optimization 

A blade structural optimization procedure 
applicable to metal and composite blades has been 
developed in which the objective function is blade 
mass with constraints on frequencies, stresses In 
the spars and In the skin, twist deformation, and 
autorotational inertia. The design variables 
(figure 8) are the total spar thickness and for 
the composite blade the percentage of ±45° plies 
(the remaining plies are assumed to be at 0°). 

This procedure and additional applications of the 
method are given in reference 10. 


5 



This section describes two example rotor 
blade designs which were developed using the 
structural design methodology. Both designs are 
based on the UH-60 Black Hawk blade. The first 
design is for a titanium single spar cross sec— 
tion. The second case has a graphite /epoxy spar 
in a single spar cross-section configuration. 

The composite spar design is compared to the metal 
spar design to explore potential weight savings 
obtained from use of the design methodology in 
conjunction with composite materials. 

Titanium cross section .- For the titanium 
spar blade , the cross - section model was based on 
the UH-60 rotor blade with identical skin, core, 
trailing edge tab, leading edge weight, and spar 
coordinates. Only the spar thickness was used as 
a design variable. The beam model representation 
of the blade used a rectangular planform similar 
to the UH-60 planform, but without any tip sweep. 

A maximum elastic torsional deformation of 3.1° is 
based on an effective aerodynamic performance con- 
straint (ref. 10). The structural constraint re- 
quires that the calculated stresses do not exceed 
the allowable material strength based on the 
Tsai-Hi 11 failure criterion. The autorotational 
capability is assumed to be the same for this 
design as it is for the UH-60. Autorotation is 
satisfied by requiring the mass moment of inertia 
to be identical to that of the UH-60 rotor system 
which is 19000 in-lbs-s per blade. Before a 
structural comparison to the UH-60 blade can be 
made, the design had to be dynamically tuned. The 
modes considered in this design were first elastic 
f lapwise and edgewise bending, first torsion, and 
second and third flapwise bending. The frequen- 
cies of these modes were required to be removed 
from integer multiples of the forcing frequency by 
0.2 per rev. 

As shown in figure 8, the minimum spar thick- 
ness needed to satisfy all the constraints was 
0.130 inch which corresponds to a blade weight of 
207 pounds. The actual UH-60 titanium spar is 
0,135 inch thick, producing a 2 10 pound blade. 

The titanium spar design is only 3 pounds differ- 
ent from the actual UH-60 blade, demonstrating 
that the mechanics of the design methodology can 
produce blade designs similar to conventional 
design processes. 

Composite cross section .- A second design was 
developed using a single T300-5208 graphite/epoxy 
D-spar. The blade models and associated design 
assumptions used in the composite design were the 
same as those used for the metal spar except for 
the spar material. Here, thickness and ply orien- 
tation of the composite spar were used as design 
variables* The plies of the spar were assuned to 
consist only of 0° and ±45* angles symmetrically 
built up. Thus, the ply orientation design vari- 
able was the percentage of ±45* plies in the lami- 
nate. The remaining plies of the laminate are 
understood to be oriented at 0°. Constraints on ' 
twist deformation, material strength, mass moment 
of inertia, and dynamic tuning were the same as 
those used for the metal design. 

Results shown in figure 8 show that the com- 
posite design satisfied the required constraints. 
Further, the minimum weight design had a 0.105 
Inch thick spar with 20 percent of the plies 
oriented at 145* degrees which resulted in blade 
weight savings of 21.5 percent. These results 


demonstrate that this design methodology, used in 
conjunction with composite materials, can result 
in significant weight savings. 


Results - Integrated Aerod ynamic Load/Dvnarai c 
Opt imitation “ 

In reference 23, an integrated aerodynamic 
load/dynamic optimization procedure was developed. 
The procedure minimized blade weight and 4/rev 
vertical hub shear for a rotor in forward flight. 
The coupling of aerodynamics and dynamics was 
accomplished by the inclusion of air load calcu- 
lations inside the optimization loop wherein the 
air loads varied with design variables. The de- 
sign model used for this procedure is the same as 
that in figure 3. The design variables include 
the stiffnesses El for spanwise and chordwise 
bending, the torsional GJ, the taper ratio, the 
root chord, radius of gyration, and nonstructural 
masses at each spanwise location. The constraints 
include upper and lower bounds on the first four 
frequencies, a lower bound on autorotational iner- 
tia, and an upper bound on centrifugal stress. 

Both single and multiple objective function formu- 
lations were used and compared. In the single 
objective function formulations, blade weight and 
4/rev shear were each individually minimized. For 
the multiple objective function formulation, a 
combination of the weight and shear was minimized 
by use of the Global Criteria Approach (ref. 11). 

A flow chart showing the logic of the optimi- 
zation procedure is shown in figure 9. The pro- 
cess is initiated by evaluating the preasaigned 
parameters (those which are constant during the 
optimization). The next step is to initialize the 
design variables and perform the blade structural 
analysis to calculate the blade properties, the 
centrifugal stress and the autorotational inertia. 
The aerodynamic and dynamic response analyses are 
performed next using CAMRAD. CAMRAD is used to 
calculate the section loads from the airfoil two- 
dimensional aerodynamic characteristics. Lifting 
line theory is used with corrections for yawed and 
three-dimensional flow effects. The blade is 
trimmed at each pass through the optimization loop 
using the wind tunnel trim option. The dynamic 
analysis In CAMRAD includes calculations of the 
frequencies and mode shapes (using a modified 
Galerkin technique) and the calculation of the 
4/rev vertical shear at the hub. A sensitivity 
analysis calculates derivatives of the objective 
function and the constraints with respect to the 
design variables. Analytical derivatives are used 
for the weight, autorotational inertia and centri- 
fugal stress. Forward finite differences are 
used for the derivatives of the hub shear and the 
frequencies. Once the sensitivity analysis is 
completed, the optimizer is called to update the 
design variables. 

The above procedure has been applied to a 
model of the Growth Black Hawk rotor blade (see 
reference 23 for details of this model). The 
baseline (analytical) model la linearly tapered 
from root to tip with a taper ratio of 3.0, has 
eight structural nodes, 14 aerodynamics segments, 
and a single airfoil for all segments. The air- 
craft is in forward flight with an advance ratio 
of 0.3. Figure 10 presents comparisons of optimum 
weight and vertical shear from the three formula- 
tions, Figure 10a compares the blade weight and 
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figure 10b compares the 4/rev vertical shear. As 
shown in the figure, the Global Criteria Approach 
provides the lightest blade structure with a sig- 
nificant hub shear reduction. This is contrary to 
the intuitive belief that the use of a multiple 
objective formulation should yield solutions lying 
between those of the single objective formula- 
tions. In other words, the blade weight obtained 
by simultaneously minimizing weight and hub shear 
might be expected to be higher than that obtained 
from weight minimization and the hub shear ob- 
tained should be higher than that obtained from 
hub shear minimization. However, this is only 
true if the objective functions are monotonically 
increasing functions of the design variables. 

This Is not true in the present case since, for 
example, the blade weight can decrease with an 
increase In taper ratio and the hub shear is a 
very complicated and nonmonotonic function of the 
design variables. 

It was of interest to determine the extent to 
which the optimization process reduced the oscil- 
latory airloads. To this end, distributions of 
vertical oscillatory airloads before and after 
optimization are compared. Figure 11 shows azi- 
muthal distributions of vertical airloads for 
the initial (reference) design and for the three 
optimum designs at a radial station of 0.75 of the 
blade span. All three optimization results indi- 
cate a significant reduction in oscillatory air- 
load amplitudes. The largest reduction occurred 
for the third formulation. As pointed out in 
reference 23, this Is a significant finding and 
indicates that the optimization process is able to 
adjust the vibration levels downward not only by 
frequency placement and modal shaping but also by 
adjusting airload distributions in a favorable 
manner. 


Concluding Remarks 

This paper has described a joint activity 
Involving NASA and Army researchers at the NASA 
Langley Research Center to develop optimization 
procedures aimed at Improving the rotor blade de- 
sign process by integrating appropriate disci- 
plines and accounting for all of the important 
interactions among the disciplines. The disci- 
plines Involved include rotor aerodynamics, rotor 
dynamics, rotor structures, airframe dynamics, and 
acoustics. The work is focused on combining the 
five key disciplines listed above in an optimiza- 
tion procedure capable of designing a rotor system 
to satisfy multidisciplinary design requirements. 

Fundamental to the plan is a three-phased 
approach. In phase 1, the disciplines of blade 
dynamics, blade aerodynamics, and blade structure 
will be closely coupled, while acoustics and air- 
frame dynamics will be decoupled and be accounted 
for as effective constraints on the design for the 
first three disciplines. In phase 2, acoustics Is 
to be integrated with the first three disciplines. 
Finally, in phase 3, airframe dynamics will be 
fully Integrated with the other four disciplines. 

This paper dealt primarily with the phase l 
approach. The paper Included: the optimization 

formulation, design variables, constraints, and 
objective function, as well as discipline interac- 
tions, analysis methods, and methods for validat- 
ing the procedure. The paper described how the 


acoustics and airframe dynamics behaviors are in- 
corporated as constraints into the design proce- 
dure. For example, acoustics imposes a local Mach 
number constraint on the blade velocity and angle 
of attack; and airframe dynamics imposes con- 
straints on the rotor blade natural frequencies 
to avoid ground resonance through coalescence of 
blade and airframe frequencies. The plan for 
validating the components of the design process 
was described and the strategy for overall valida- 
tion of the design methodology was defined. These 
validations are critical to the success of the 
activity and are viewed as the primary products of 
the work. Finally, some representative results 
from work performed to date are shown. These 
include aerodynamic optimization results for 
performance, optimal placement of tuning mass for 
reduction of blade shear forces, blade structural 
optimization for weight minimization subject to 
strength constraints, and integrated airload/ 
dynamic optimization results for vibration 
reduction. 

The results of the individual optimization 
procedures demonstrate the potential of optimiza- 
tion in design of future rotorcraft, both from the 
standpoint of efficiency of the process as well as 
potentially improved products. The results of the 
integrated airload/dynamic optimization procedure 
demonstrates that there are significant benefits 
awaiting analytical designers who pursue interdis- 
ciplinary design approaches. 
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TABLE 1 SUMMARY OF DESIGN VARIABLES 


Description Symbol 


Tuning mass at location 1 

Spanwlse location of i-th mass x^ 

Wing box dimensions tj, t2» t3 

Ply thicknesses C$5, c 0 

Depth of blade at root h r 

Ratio of blade depths at tip and root - h r /h c 

Maximum pre-twist of blade T max 

Percent blade span where taper begins r 

Blade root chord c r 

Airfoil distribution 

Hinge offset e 

Blade angular velocity fl 

Number of blades on rotor N 

Blade radius r 

Ratio of root chord to tip chord X c » c r /c t 
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TABLE 2 SUMMARY OF CONSTRAINTS 


Constraint Description 

Form of Constraint 

Comments 

Main rotor horsepower 

HP i < HP avail for 

For 5 flight 


i-th condition 

conditions 

Airfoil section stall 

C D * C Dmax 

Enforced at 
12 azimuthal 
locat Ions 

Blade frequencies 

f . , < f , < f, 

^ < 1 v lu 

v ik max 


Blade vertical load 


Blade inplane load 

H ik <>.x 


Transmitted in-plane 

X. < X 

Y k < Y maX 
T k max 


hub shears 


Hub pitching moment 

P k < P max 


Hub rolling moment 

P k * ^max 


Blade response amp. 

^k * ^ax 


Autorotat ional inertia 

Eaii r i > a 


Aeroelastic stability 

Re (\) < -e 


Wing box stresses 

R < 1 

R - Tsai-Hil 1 
criterion 

Blade tip deflection 

w < w 

e < e 

max 


Blade twist 


Blade tip Mach no. 

M < M m „x 
h < V.X 

Limits 

Blade thickness 

thickness 

noise 

Blade lift distribution 

dC l /dx < S «ax 

Limits BVI 
& loading 
noise 

Ground resonance 

l Q ' “ill < “af 
*1 < « 4 f u 

Effect ive 

Rotor/Airf rame 

airframe 

frequency coupling 

constraints 


TABLE 3 INTERACTIONS AMONG DISCIPLINES 


Variable 

Acoustics 

Aerodyn. 
(Perf & Loads) 

Dynamics 

Structures 

Fuselage 

Dynamics 

Airfoil Dist. 

S 

S 

W 

W 

W 

Planform 

S 

s 

S 

S 

s/w 

Twist 

W 

s 

S 

W 

w 

Tip speed 

s 

s 

S 

S 

s 

Blade number 

s 

w 

S 

w 

s 

Stiffness 

w 

s 

s 

s 

s/w 

Mass dist* 

w 

w 

s 

s 

s/w 

Hinge offset 

w 

w 

s/w 

w 

s/w 


S - Strong interaction 
W - Weak interaction 
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TABLE 4 CANDIDATE TASK AND MISSION 
FOR PHASE 1 DESIGN ACTIVITY 


Description 

Specification 

Condition 

4000 ft 95* 

Aircraft gross weight 

16875 lb 

Installed power limit 

3400 HP 

v crulse 

140 kts 

V 

v max 

200 kts 

g’s at 120 kts 

3.5 

Vertical rate of climb 

1000 fpm 

Airframe structure 

UH-60B 


Other constraints and guidelines are 
specified in table 2. 


Phase 1 

Optimizer i 


| Aerodynamics 
, Blade dynamics 
! Structures 

Acoustics 
AMrevne dyn 

Constraints 

only 


j Phase 2 

J Optimizer 

Aerodynamics 
Blade dynamics 
Structures 
Acoustics 

Airframe dyn 

Constraints 

only 


Phase 3 

V Optimizer — 

! Aerodynamics 
j Blade dynamics 
| Structures 
| Acoustics 
j Airframe dynamics 


Fig. 1 Phased approach to development of 
Integrated rotorcraft optimization 
procedures . 



Fig. 2 Integrated rotorcraft optimization 
development plan. 




Current dssign variables L 
[ and fixtd paramstara 


Aarodyrn 

anafys 


(load*) 


i 


Paaign vartabta praprocattoft 


Airload# h 


Dynamic 

analyiia 


Dynamic] Structural 
load* analysts 


Rasponaa 

stability 


Aerodynamic j j ( 

anafyafa —►Objactive function and constraints u 
(parforman ca) [ ^ ’ 

HP rag'd 


Sensitivity analysis 


; Strassas 

'deflections 


Optimizer 


i Updated design variable* 


Fig. 4 Integrated aerodynamic-dynamic-structural 
optimization of rotor blades. 





Objective 


Design 

variabias 


Dasign lima 


I 

l 

1558 hp 

1533 hp 

Twist, deg 

-12 

-15 

Percent taper 

.80 

.91 

Taper ratio 

3.0 

3.1 

Root chord, ft 

2.3 

1.78 


5 weeks 

2 days 


Fig. 5 Results of aerodynamic performance 
optimization. 
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• Design goal * Find optimum combination of 
masses and their locations to reduce blade 
root vertical shear 

• Method * Formulate optimization procedure 

* Use masses and locations as design variables 

• Minimize - 

• Blade root vertical shear 

• Added mass 


Fig. 6 Selection of optimum locations of tuning 
masses for vibration reduction. 



Fig. 7 Time history of vertical root shear 
minimized for 2 modes/3 harmonics. 
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Twist constraint 

8«.i 


Autorotation constraints 



Material strength constraints 



R = Tsai-Hill 
failure criterion 


Results 

! Parameter Actual Metal Comp.} 

i UH-60A spar s par 

Spar mat. Ti Ti GnE 

.Spar thick., in. 0.135 0.130 0.105 

45'plys, % — — 20 

Weight, lbs.. 210 207 163 

Margin (1-R) — 0.103 0.000 

I Twist, de g. ™ 0.92 2 55 


% Of 45° 
plys (all 
others 0 s ) 



Design variables 
t (thickness 
of spar) 


Fig. 8 Structural optimization for minimum weight 
rotor blades. 


Hub shear 
(Ib/fl) 



(b) Optimum 4/rev hub shear from various 
f ormulat ions . 



j 


Fig. 9 Flow chart for integrated aerodynamic 
load/dynamic optimization procedure. 


Fig. 10 Concluded. 


o Reference 
a Weight 



Fig. 11 Azimuthal distributions of vertical 
airload, 75Z radius, p * 0.3. 


Weight 

(lbs) 



(a) Optimum weight from various formulations. 

Fig. 10 Results from integrated aerodynamic load/ 
dynamic optimization procedure. 
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